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Abstract 

In Split-Supersymmetry models, where the only non-Standard Model states 
produceable at LHC-energies consist of a gluino plus neutralinos and charginos, 
it is conventionally accepted that only mass differences among these latter 
are measureable at the LHC. The present work shows that application of a 
simple 'Kinematic Selection' technique allows full reconstruction of neutralino 
and chargino masses from one event, in principle. A Monte Carlo simulation 
demonstrates the feasibilty of using this technique at the LHC. 
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1 Introduction 



As data from the LHC is recorded and analyzed in the upcoming years, experi- 
mentalists will look for signatures of physics beyond the Standard Model (SM), in 
particular Supersymmetry (SUSY). The most well-studied scenario, in which SUSY 
alleviates the hierarchy problem with a low-energy (sub-TeV) spectrum of sparticles, 
entails copious production of strongly-interacting squarks and gluinos, identified via 
their associated jets [H[2], which cascade through numerous decay channels involving 
other sparticles, including e.g. electroweak (EW)-interacting sleptons, neutralinos, 
and charginos — the masses of these sparticles, whose precise values are crucial to 
understanding features of an underlying fundamental theory, may be reconstructed 
(at least partially) from measurements of various invariant mass endpoints in certain 
exclusive decay channels (e.g. [3]). Of course it is also entirely possible that the SUSY 
spectrum is far above the TeV level, hence inaccessible at the LHC. In between the 
above two extremes, phenomenologically speaking, is the scenario where some of the 
sparticles are light, while others extremely massive and decoupled, "Split SUSY" [I] 
providing the most popular example. 

At low energies, the Split-SUSY spectrum, aside from the established SM particles 
and one light Higgs boson, contains only four neutralinos i = 1..4), two charginos 
(Xj > j = 1) 2), and a long-lived gluino (g). The phenomenology of this latter, which 
would be expected to form so-called "R-hadrons" in the detector (and therefore not 
immediately decaying to other sparticles) has been thoroughly covered elsewhere 
(see [E1E]) and will not concern us her^|. The focus of the present study is rather on 
the neutralinos and charginos, hereafter collectively referred to as 'EW-inos'. These 
cannot decay via squarks and sleptons, which are many orders of magnitude heavier, 
but must rather decay promptly via a Higgs or EW gauge bosons (Z°, W^). If 
mass differences between EW-inos are smaller than mz or rriw, they will undergo 
3-body decays to quarks or leptons: e.g. xi ~^ Z°*(— > ^i^, qq)Xii xt ~^ Z *(^> 
qq)Xij an d Xi ~> W ± *(—> ^v, qq')x.i- I n particular, dilepton pairs from 
the above off-shell Z decays will have an invariant mass distribution which cuts off 
sharply at m^o — m^o or m~± — m~± , so the usual conclusion from Split-SUSY studies 
is that only mass differences between EW-inos are measurable at the LHC [6]. 

The thesis of this work is that one can do much better than just find EW-ino mass 
differences — the masses themselves can be reconstructed from additional kinematic 
analysis made possible by the fact that EW-inos must be pair-produced to preserve 
R-parity, which is theoretically motivated in making the lightest SUSY particle (LSP) 
a good dark matter candidate. Each SUSY EW-ino event (x^Xj^ > Xi x/ \ or Xi°Xj°) 
may thus contain multiple (as many as ten) hard leptons, the momenta of which, 
when contracted into all possible invariant masses (as in [7]), encode much informa- 
tion. This was previously overlooked in the literature, presumably because pairs of 
EW-inos arising from hadronic collisions carry an uncertain center-of-mass energy, 
hence yielding final state leptons not amenable to the usual invariant-mass endpoint 

2 Even if the gluino does decay in some corner of parameter space, this will only assist with the 
current study by boosting signal rates. 



analysis. As shown in Hidden Threshold (HT) methods [8], however, correlations 
among such invariants (i.e. a Dalitz plot) still carry information about endpoints 
and, more importantly for the present work, distribute events according to the kine- 
matics in each respective EW-ino decay frame. The strategy in this work, therefore, 
is to focus on one region of a Dalitz plot where events must arise from the same 
decay frame kinematics, find the Lorentz-boosts back to the frames of the decaying 
EW-inos, and match energies/momenta (including measured missing transverse mo- 
menta) to extract the values of relevant masses. A single (perfect) event may suffice 
for full reconstruction in principle, though in practice (including detector effects and 
backgrounds) one must do with a collection of less-than-perfect events which will 
give statistical distributions of the unknown masses. 

In the following, let us then proceed thusly: Section [2] will explain this 'Kinematic 
Selection' method in the context of Split-SUSY EW-ino decays; Section [3] will detail 
how to reconstruct EW-ino masses from a perfect event and Section H] will then test 
and confirm the feasibility of this in a Monte Carlo simulation appropriate to the 
LHC environment. Summary and comments on further elaborations are contained 
in Section O 



2 Kinematic Selection Technique for EW-inos 

In Split-SUSY models, EW-inos at the LHC can only be pair-produced in quark- 
quark s-channel processes through an off-shell W, Z, or 7: 

qq' — > W* - xt^ (1) 
qq^Z* - x-X? (2) 
qq^Z*/j* - xtxj (3) 

Then, since each EW-ino cannot decay through squarks or sleptons, but only through 
a Z° or W ± (or a light higgs h°, though this tends to be subdominant), it must decay 
among the following five tree-leveo channels (taking m~± > m~o and m~± > m~o): 
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where the Z° or W ± could be on- or off-shell. The number of possible decay chains 
combining ([I])-© and (jl])-©, even without distinguishing on- or off-shell interme- 
diaries or considering the rest of the decay chain, is already quite large, but most of 
these, fortunately, will not be needed in the present study. 

3 Loop- level decays can also have phenomenological importance, e.g. Xi° — > 7x1° [9]. 



2.1 Chargino-Neutralino Modes 



The most heavily-produced state in Split SUSY models is likely to be a chargino- 
neutralino pair xtxli srnce these sparticles are relatively light and well-mixed, where 
it is further assumed that they proceed to decay through an offshelE Z° or W ± to 
leptons (£ = e, /z): 

xt - W^i^tv)^ (9) 

xl - z°*H^ T )x? (io) 

The endstate will therefore contain three leptons (of which at least two are opposite- 
sign-same-flavor (OSSF)) whose momenta ^1,2,3 can, in the spirit of [7], be system- 
atically contracted into three independent invariant masses^]: 

M\ EE (P1+P2+P3? (11) 
EE {(pi+p 2 -p 3 ) 4 +(P3+Pl-P2) 4 + (P2+P3-Pl) 4 }/3 (12) 

M\ ee {( Pl+P2 y + ( Pl+ p^ + (p 2+ p 3 y}/3 (13) 

The problem now, in which the HT technique assists, is how to use the infor- 
mation contained in the above invariants to select events with a desired kinematic 
configuration. Observe, first of all, that the off-shell W in (pQ) will itself have an in- 
variant mass somewhere (and unpredictably) in the range m~± + m^o < m* < Elhc 
where E^hc ~ 14TeV is the theoretical maximum partonic collision energy at the 
LHC. Let us consider the case where m* = m~± + m^o and designate this 'threshold 

Xl A2 ° 

production'. 

To simplify the discussion, assume we have a e + e~fi ± endstate (the following 
will also pertain to same flavor states e + e~e ± or fj, + fi~fj, ± with the correct lepton- 
pairing). From relativistic kinematics, it is quite straigtforward to show that, for 
threshold production, when M e + e - = {p e + +p e -) 2 is maximal, M\% is minimal when 
the kinematical configuration in Fig. [Th. is attained: in the rest frame of the par- 
ents xt and X21 the electron and positron are produced back-to-back with maximal 
momentum along directions perpendicular to the muon, which also carries maximal 
momentum (see Appendix for derivation). The minimal value of M^z is then given 
by 
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(the other invariants My and M« are, on the contrary, trivially minimized by p^± = 0, 

giving M™ in = M e T e - and M™" = (l/3) 1/4 M e T e - , respectively). For the more 
realistic case of non-threshold production, i.e. m* > fn~± + m^o giving a relative 



4 If decays occur through on-shell Z° and M /± , the signal is much more challenging to extract, 
being swamped by WZ and ZZ backgrounds. 

5 These have the advantage of systematic definition and symmetry under lepton interchange at 
the cost of algebraic complexity. 



velocity (3 between the Xi and X.2i ^ i s a l so easy to show that M e + e - is maximal 
and M121 minimal (for this (3) for the same decay configuration (in the respective 
xf and X2 'parent-frames') of Fig. The relevant Dalitz plot is therefore "M e + e - 
vs. M121" , the events of interest accumulating along where the line M e + e - = M™ ^ 
(= m^o — m^o) hits the kinematically-allowed portion of the ploi@. 




Figure 1: Schematic of the kinematic configurations (with e «-> /i as well) which mini- 
mize various invariant masses at threshold for (a) Chargino-Neutralino and (b) Neutralino- 
Neutralino modes (these particles decay at rest in this frame). 



2.2 Neutralino-Neutralino Modes 

Neutralino pair production (j2J) may also be significant, but only for unlike neutralinos, 
i.e. i 7^ j, due to a suppression of the Z°x^X4 coupling [HI [TO]- Assuming again that 
the neutralinos decay through a 3-body decay as in (ITU1) . we now have an endstate 
described by four lepton momenta ^1,2,3,4 which can be analyzed via a set of seven 
invariant masses [7], including e.g. 

Ml = (p! + p 2 + p 3 + p A ) 2 (15) 

M 2l2i = {(Pl+P2-P3-P4) 4 +(Pl+P4-P2-Ps) 4 + (P2+P4-Pl-P3) 4 }/3 

Going through the same argument above for X1X2 modes, one finds that the threshold 
kinematic configuration in Fig. [lb, where M e + e - is maximal, forces M4; and M 2 m to 
attain the minima (derived in Appendix) 
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yj(rng) — m^o)(m^> + m^p — m%a — m-o/m^q) (16) 
2(m^o — m^o) 2 + (m^o — m^o + m^o — m~ / m^o) 2 ^ 

v X-i X\' v Xj Xi Ai * Xj ' 



and events of this type can be found near these minima on a plot of M e + e - (or M^+^- ) 
versus M Ai and M 2 ni- 



6 One might hope to measure the endpoint p4[) from the intersection and thereby constrain SUSY 
masses, but this turns out to require very high (sub-GeV) endpoint precision and is much inferior 
to the present method. 



Figure 2: Class of parent- frame kinematics sought: in the frame where the x\ decays at 
rest, the e 1 * 1 are along global coordinates ±y with maximal energy; likewise in the X\ decay 
frame the muon is along ±x with maximal energy. Any velocities f5\^ of these decaying 
states in the lab frame are permitted. 



3 Finding Masses from One Event 

Let us now suppose we have a e + e~/i ± event of maximal M e + e - and minimal M^z 
(a very similar analysis can be done for e + e~fi + fi~ events from Xi X°j decays, but we 
will not need this for the present stud}0). By the discussion of the last section the 
kinematic configuration must be of the type shown in Fig. [21 i.e. if the decaying x\ 
and Xi na d no motion then the would have equal and opposite momenta (along 
say ±y) while the muon would be emitted perpendicular to y, say along x, also with 
maximal kinetic energy; the X2 an d xf are ; however, permitted to be moving with 
different velocities 3 1,2, so the observed leptonic momenta will generally point in 
random directions. 

These leptonic momenta nevertheless carry useful information, for if we knew 
3 1 9 as well we could reconstruct all three unknown masses m^o , m^o , and m~± as 

follows: the observed missing transverse momentum ^>t must clearly arise from the 
invisible particles 2xi + ^, and since their 4-momenta are known functions of the 
masses in the X2 anc ^ Xi decay frames, it is a simple matter of Lorentz-boosting 
these by (3 \^ to the lab frame and matching to the two components of fix- Two 
matching conditions plus the dilepton edge (M^L = m^o —m^o) determines the set 
{rn^o, m^o, m~±}. 

It is easy, in fact, to determine 3 1: this corresponds to the unique Lorentz 
transformation Ai which makes the transformed e momenta equal and opposite 
(= ±p'), as well as simultaneously bringing the corresponding Xi to rest, and is 
given by 

% = + I*' (17) 

E e + + E e - 

As for 3 2, there is no nice analytical expression, but we can nevertheless constrain 



7 If the number of e + e n + u events is sufficiently high we can also get events where M e + e - and 
are both maximal, also allowing for a straightforward mass-reconstruction 



it by conservation of the total missing 4-momentum 
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(ie. total missing 4-momentum = 4-momentum of one LSP plus 4-momentum of 
other LSP+i/ system) of which the two transverse components fa are measurable, in 



addition to the three kinematic constraints 
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Thus, ( 1151) and f[T9l compose a system of five equations for the six unknowns { (3 2, 
m^o, m^o, m~±}. If, as in a wide class of SUSY spectra (here controlled by SUSY 
parameters fi, M± t 2, and tan(3) we have the approximate relation m~± « m^o, then 
we are less one unknown and the five equations can be numerically solved for the 
masses m 7 o and m 7 o. 



4 Monte Carlo Test 

In a real experiment, there are of course many reasons why even a 'perfect' event like 
Fig. [2] does not suffice for reliable mass reconstruction: measurement errors as well 
as inherent finiteness of detector resolution and sparticle widths will throw off the 
solution. Then there is the reality that no event is perfectly situated at an endpoint 
and, moreover, competition from backgrounds is expected. What we must do in 
practice, therefore, is to collect a number of events in some optimized neighborhood 
of the region of interest on the "M^- vs. M; 2 /" plot, impose conditions (fITl) - (Tl9l) 
on each event (with M^L — ► Mt+i- plus the assumption m~± w m^o ), and study 
the distribution of extracted masses m^o 2 (the reader is referred to the Appendix for 
a more detailed discussion of the numerical solving procedure). 

Let us see how this might work by running the above programme through a Monte 
Carlo simulation of LHC data. Suppose for defmiteness that Nature has chosen the 
Split SUSY parameter point considered in [6] with GUT-scale parameters 

M 1 = M 2 = M 3 = 120 GeV 
= -90 GeV 

tanf3 = 4 

in addition to a symmetry-breaking scale m = 10 9 GeV. Integrating down to EW 
energies (Q = mz), all SUSY particles decouple except for the gluino and EW-inos, 
which attain the spectrum shown in Table 1; this is consistent with LEP and dark- 
matter constraints. At LHC energies the dominant chargino-neutralino production 
channels would then be^l xtx.2 ( a = 4650 fb) and xfxt ( a — 2099 fb), while the 
main neutralino-neutralino channel is X2X3 { a = 876 fb). 

8 See [5] for a complete table of these. 



Table 1: Relevant masses (in GeV) at the Split- SUSY point we consider. 
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LHC EW-ino events (pp -> xt 2 X°, Xi^Xi^) and SM backgrounds ZZ, IfZ 

and Wj* (see p2] and [33] for a good discussion of these and others not necessary 
for this study), corresponding to 300 fb~ l integrated luminosity are then generated 
via the HERWIG 6.5 package [H] and run through a simplified detector simulatoi0. 
The following cuts are then employed, depending on the number of final leptons: 

For 2-Lepton Endstates: 

• Leptons must be isolated: no tracks of other charged particles are present in a 
r = 0.3 rad cone around the lepton, with less than 3 GeV of energy deposited 
into the electromagnetic calorimeter for 0.05 rad < r < 0.3 rad around the 
lepton. 

• Leptons must be sufficiently hard: p l T > 10, 8 GeV for £ = e, \i. 

For 3- and 4- Lepton Endstates: 

• Leptons must be hard and isolated as for 2-lepton endstates. 

• Sufficient missing energy must be present in each event: Efr > 50 GeV. 

(note: J^-lepton selection criteria are shown for completeness only; the present study 
does not consider them further due to smallness of rate at this Split SUSY point) 

Though SM backgrounds are substantially reduced by these cuts, they still tend 
to far outnumber the SUSY signal events for 2- and 3-lepton endstates. We will 
see shortly there is no cause for worry, however, since SM backgrounds populate 
uninteresting regions of the relevant invariant mass plots, shown in Fig. [3l 

First observe the large number of 2-lepton events in Fig. [3^: there are roughly 
6 ■ 10 4 SUSY £ + £~ events plus 3 ■ 10 5 SM events (mostly from Z decays and hence 
sitting near the Z-pole, Mp+p- ~ 91 ± 10 GeV, not shown in the Figure) which, after 
subtracting wrong-flavor e ± fi ± combinations (5 • 10 4 of these total), give us a dilepton 
invariant mass distribution that clearly identifies an endpoint at Mp+p- ~ 39 GeV due 
t° X.2 ~^ X? decays^ (a second endpoint from Xi decays near ~ 70 GeV is barely 

9 The set-up is the same as in several of the author's previous publications (e.g. [71 [TUl [23 EE] ) , an d 
includes a privately-coded fast detector response simulation incorporating all the requisite simplified- 
geometry calorimetry, missing energy reconstruction, lepton isolation, etc., and has been checked 
against results in the literature using publicly available codes. 

10 In this feasibility-level study, it is sufficient to mark this endpoint to within a few GeV; more 
complete analyses [TJ [2] of dilepton distributions with comparable or lower statistics verify harm- 
lessness of SM backgrounds and suggest sub-GeV level precision is easily attainable. 




Figure 3: Invariant mass plots for 300 fb^ 1 luminosity (SUSY + SM): (a) The flavor- 
subtracted dilepton invariant mass distribution clearly identifies the endpoint at ~ 39 GeV. 
(b) correlation between M e +£- and M^z; events are taken in the boxed region shown, i.e. 
the neighborhood of where the line M^+£- = 39 GeV hits the envelope. 

discernible but might be claimed statistically significant via a more in-depth analysis, 
e.g. 17)0 

Turning now to 3-lepton events, SUSY events (e + e~/i ± + /i + /i~e ± as well as same- 
flavor e + e~e + /i + fj,~ /r* 1 event 9^1) number close to ~ 8000 against a SM background 
several times larger, but this latter is, in the 'Mf+f- vs. plot shown in Fig. [3b, 

concentrated mostly up near the Z-pole again, and to a much lesser extent throughout 
the bulk of the plot. The kinematically allowed region has a fairly clear diagonal 'left 
edge', and it is where the line M £ + £ - = 39 GeV hits this edge that we should expect 
events of the type in Fig. [2j Events are therefore collected from the boxed region 
shown, the limits of which (M^- = 39 ± 4 GeV, M; 2 ; < 40 GeV) give the optimal 
distribution of m^o, shown in Fig. UK Although this region certainly includes a large 
number of background events (e.g. X2X3 ~~ > e+e ~ H + H ^Xii with one of the four 
leptons failing the hardness cut, or X2 X^t - * e+e ~/ i=Fz/ Xi)) SM processes as 
well as same-flavor signal events with the wrong lepton-pairing) , their presence can 
be tolerated since these either do not give a physically-acceptable solution to the 
system of equations f|T7l) - ffT9l . hence are rejected, or they give no preferred solution 
and lead to a uniform 'noise' in the m^o-distribution. This latter, in fact, peaks 
sharply (a Gaussian fit gives m^o ~ 63 ± 3 GeV) but somewhat lower (by about 
15%) than the nominal value in Table 1 (the m^o-distribution looks the same but 
shifted by ~ 39 GeV). This deviation may arise from the fact that the assumption 
m~± = m^o is inaccurate by several percent at this parameter point (see Appendix for 



11 Thcre is also an 'edge' at Mp+i- ~ 20 GeV, but this is actually an effect of lcpton pT cuts: 
since two leptons with momenta p± and relative angle 6 give an invariant mass of Mi+ e - = 
^/2p + p_(l — cos0), cutting away p± < 10 GeV will tend to deplete the M t + t - spectrum below 
2p ± = 20 GeV. 

12 For same-flavor events, both possible lepton-pairings are plotted. 



a discussion). If, instead, one does not make this assumption, but knows beforehand 
the value of m^o from other measurements, (1TTI) - (1T91 can also be used to find ni~±. 
The distribution shown in Fig. Sb gives the correct m~± within errors (tu^± = 108 ± 
15GeV). Note this value of m~± can be put back into (|T7|) - (|T9|) to solve for m^o 
again, iterating the process. 
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Figure 4: (a) Distribution of m^o from events in the boxed region of Fig.\^j), assuming 
m~± = m~o. (b) Distribution of m~± assuming the correct value of m~a = 71 GeV. In both 
plots the LEP2 constraint m~± > 100 GeV fTSlj has been applied. 



5 Discussion and Conclusion 

This paper has introduced a Kinematic Selection technique applicable to EW-ino 
pair decays where 3- or 4-lepton events with specific parent-frame kinematics are 
captured for analysis. Such a technique is particularly suited to Split SUSY models, 
where EW-inos are the only low-lying states expected to be observable at the LHC. 
At the specific parameter point studied, the rate of xtx2'P a ^ T decays to three leptons 
was sufficient to capture O(10 2 ) interesting events, the kinematics of which could be 
reconstructed well enough to extract two EW-ino masses m T o within 15% or so of 
their nominal values, assuming that m~± = m^o. It is, of course, entirely possible 
that other Split-SUSY parameter points would give higher rates and thus allow us 
to study X2X3 4 modes as well as those involving the heavier chargino xf (here the 
kinematics are more complicated, but the final state also contains more leptons, hence 
more useful invariant mass constraints). Also, if EW-inos decay through an on-shell Z 
or W, or through a light Higgs boson, we can in principle look for leptonic invariant 
mass correlations which isolate events with specific parent-frame kinematics (note 
final state jets can also be included in this formalism); the method is quite flexible. 

What other mass reconstruction methods are available for analyzing Split-SUSY 
EW-ino decays? First consider the neutralinos. There are now an array of methods 



which take advantage of the pair-production of neutralinos. One class of "Mass-Shell 
Techniques" (MST), represented in the work of pLQ] and [2D], essentially depends 
on maximizing the solvability of assumed mass-shell constraints in a given sam- 
ple of events. This seems quite effective for on-shell decays^!, but for the off-shell 
decay topologies in the present work these methods cannot be applied since there 
are not enough such constraints. Recently fashionable "transverse mass variable" 
methods [221 [23], e.g. m T2 , might be applied, though these are usually stated for 
symmetric decays. In one such development [24J, for example, a "constrained mass 
variable" m 2 c proves quite powerful for X2X2 modes (followed by off-shell decays such 
as ffTUj) ); though such modes are expected to be negligible in Split-SUSY scenarios, 
presumably m 2 c could be applied to the case of unlike neutralinos X2X3 as well- 
It's worth mentioning here that in the some of the latest developments with m^, 
e.g. mr2-Assisted-On-Shell (MAOS) reconstruction [25], information on the full LSP 
4-momentum can be gleaned for both mass and spin determination. 

As for decay modes with charginos such as X1X21 author is not aware of any 
work showing how to reconstruct all the unknown masses (MAOS has not yet been 
tested [2S]) — perhaps the above techniques can encompass these modes as well, but 
there may be fundamental difficulties with extra invisible particles (neutrinos) in the 
decay products (e.g. MSTs would have too many unknown degrees of freedom in each 
event). Finally, there is the under- addressed question of multiple competing decay 
channels, e.g. when several different XiX.'j an d XiX°j occur with similar rates. The 
case of XiXj yielding a e + e~fi + fi~ endstate, in particular, is subject to a wedgebox 
analysis [15] to partially separate events according to decay topology, though this has 
not been (but should be) extensively tested for mass reconstruction methods which 
have so far only concentrated on a single channel. Note that in the method of the 
current paper this separation is unnecessary, since minima such as M^™ for various 
channels lie on different points of the envelope in Fig. [3b. It seems quite natural that 
a combination of several techniques will be necessary to both isolate relatively pure 
samples of a given decay and reconstruct unknown masses as best as can be done at 
the LHC For example, an MST-analysis might be applied to 4-lepton events from 
X2X3 modes to get a ballpark estimate of m^o, this value then used in the current 
method with xtx.2 modes to determine rn^±, as done in Fig. Sb. 

In conclusion, then, this work represents the first application of a Kinematic Se- 
lection technique, found to be of particular use in Split-SUSY models. The strengths 
of Kinematic Selection include simplicity (relativistic kinematics) and robustness 
(works for multiple decay channels, even with backgrounds), which should make it a 
useful tool to experimentalists unraveling data from the LHC. 
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But see [21] for some important caveats. 



Appendix 



Derivation of Minima 



In Section [2] the threshold minima (|14p and ffT6l) were stated without derivation; here 
let us see how these were obtained. 

Starting with M121 for a e + e~/i endstate, choose the threshold frame of reference 
to be such that the electron/positron (of maximal invariant mass) are emitted along 
the z-axis ±z, while the muon (with maximal energy) is produced at spherical angles 
0. 6: 
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Plugging these four- vectors into the definition of M121 in ( {TBI and simplifying a bit, 
one obtains 
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This is clearly minimal when 9 = 7r/2, and further algebraic simplification leads to 

With a four-lepton endstate like e + e , where the e + e pair, say, has maxi- 

mal invariant mass along the z-axis, one muon (fi + ) will be going at an angle 9 + to 
the z-axis, while the other (/z~) has its own spherical angles 6L, (f>: 
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where the energies E± take on a range of values determined by the relative angle 
between the muons and 3-body kinematics. One could then plug these expressions 
into the definitions of, e.g., M41 and M 2 m and minimize over the angles, but it's much 
faster to intuit that since we're interested in minimizing an invariant mass, the muons 
should be going in the same direction (9 + = 9_ and = 0) which forces E + = E_ = 
(m~ — m- )/2myo; plugging this into the definitions ffTBT) and simplifying yields the 
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quoted minima f|T6|) . and numerically sampling over 
are indeed correct. 



#_,0)-space confirms these 



Numerically Solving for the Masses 

The most direct method of numerically solving the system of equations (lisp and (1191) 
for the variables {/?2x, ^2y, ftiz-, m x?> m x2' m xt ^ ( wnere > sa y> m x? * s known ) is to 



simply loop over a liberal range of values for m~± (say between 50 GeV and 300 GeV 

) and the components of j3 2 (02x,y,z e&ch between -1 and 1, with | (3 2 \ < 1), imposing 
the other constraints inside these four nested loops. In fact, we only need to loop 
over two of the components of /3 2 , since the third is fixed by a requirement on the 
transformed muon energy: when the muon four-momentum p v is boosted back to the 
chargino's decay frame by /3 2 , i.e. 
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we must satisfy the constraint from 3-body kinematics, 

E - P 2 xPx ~ P2yPy 
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This can be rearranged into a quadratic equation for j3 2x , 
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(20) 



There is of course a potentially two-fold ambiguity in the solution for f3 2x , and both 
values must be tried (if they are indeed real and satisfy \f3 2 \ < 1). Scanning over the 
3-dimensional (m~±, /?2j/, /?2 Z ) -space then, we look for the point which best satisfies 
the missing-momentum constraints (these must be satisfied within ±10 GeV) as well 
as the kinematic constraint p ' ■ p " — 0, achieved by minimizing A: 

A ^ y/A* + A§ + Af 

Ai = f + Pixim^o + 

P2x72(m~ f - E") - (1 + ( 72 - \)%p" x - (72 

P2 

A 2 = f + Pxy^m^o + 

[3 2yl2 (m~ t - E") - (1 + (72 - - (72 

A 3 = a(^'-^")/(\^'\\^"\) 

where (3\ is already known from fllTj) . and a is a weight (high a ~ 1000 seems best, 
meaning that all solutions have essentially perpendicular leptons, | cos^l < 0.01): this 
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minimization always gave a unique solution in all cases tested. Moreover, this pro- 
cedure is efficient in dealing with backgrounds (or same-flavor e + e~e ± and /i + /i _ /i ± 
events with the wrong lepton-pairing) : either these fail to satisfy both missing energy 
constraints (Ai^) within ±10 GeV, or yield solutions randomly distributed across 
mass space, which merely provides a uniform 'noise' in the solution histogram. 

The same algorithm applies, of course, when m^o is not known (so we loop over 
it) but we assume m~± = m^o. Here, however, since actual kinematic data comes 

' Xl A2 

from an event where this equality does not strictly hold, we expect the result of 
the numerical solution above to have a systematic error: setting m~± = m^o + e, 
the numerical solution m'~ Q is offset from the actual value, m'~ = m^o + 5. To 

Xi ' Xi Xi 1 

quantitatively understand the relationship between 5 and e, we would be best off 
running many simulations at different Split-SUSY parameter points and plotting the 
correlation '5-versus-e'; but since this is extremely time-intensive and not practical 
for the present work, we can get a quick-and-rough idea in the perturbative limit 
(5, e « mzo) by letting (m'~ , mL±, m'~ ) and (m^o, m~±, m^o) both solve Ai = 
for the same Lorentz boosts, where 

m~± = m^o + e (21) 
m'L = ml + 6 
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and then seeing what the relationship between 5 and e must be. Thus, the exact and 
approximate solutions, respectively, satisfy 
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Explicitly inserting the 5- and e- dependencies from (1211) into (1221) and ( 1231) and setting 
these latter equal, 
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And now expanding and keeping only terms of 0(5) or 0(e), we finally arrive at 
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(24) 



From fl24l) we see that, even in this greatly simplified treatment, the sign and magni- 
tude of the correlation between 5 and e depends on the Lorentz boosts (3 1,2 specific 
to each event. All things being equal, however, fixxPfx ~ —$2x12 (the x° an d xf ma J 
tend to go in opposite directions), we can drop the smallish r-dependent terms, and 
we might thus expect a negative correlation 5 ~ — e-. This is indeed what is observed 
at the Split-SUSY point in the present study, where S ~ — 8 GeV and e ~ 5 GeV. 
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